Cold Nuclear Matter Effects on J/tp Yields as a Function of Rapidity and Nuclear 
Geometry in d+A Collisions at v /sJ~^~=200 GeV 
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We present measurements of J/tj) yields in d+Au collisions at ^/s NN =200 GeV recorded by the 
PHENIX experiment and compare with yields in p + p collisions at the same energy per nucleon- 
nucleon collision. The measurements cover a large kinematic range in J/ip rapidity (—2.2 < y < 2.4) 
with high statistical precision and are compared with two theoretical models: one with nuclear 
shadowing combined with final state breakup and one with coherent gluon saturation effects. In 
order to remove model dependent systematic uncertainties we also compare the data to a simple 
geometric model. We find that calculations where the nuclear modification is linear or exponential 
in the density weighted longitudinal thickness are difficult to reconcile with the forward rapidity 
data. 



PACS numbers: 25.75.Dw 

The measured yields of quarkonia states in p+A (or 
d+A) collisions provide information about the time scale 
and dynamics for the creation of a cc pair and its evo- 
lution to a color-singlet quarkonium state. The prop- 
agation time of the cc pair through the nucleus is set 
by the incident energy of the projectile and target and 
by the relative longitudinal momentum of the cc pair. 
Fixed target p+A experiments at Fermilab [l[ reveal a 
substantial suppression for forward rapidity J /if) and ip' 
at a similar level, leading to the conclusion that the sup- 
pression must occur at the prehadronic stage. An anal- 
ysis of results for -Js NN — 17-42 GeV indicates that 
in addition to modified initial production due to nuclear- 
modified parton distribution functions (nPDFs) , a break 
up cross section (ov) for the cc precursor state to the 
J/ip is important, and that o^r decreases as the relative 
center-of-mass energy between the cc and the nucleon in- 
creases. Extending these results to collider energies at 
RHIC is important. The dominant production mecha- 
nism for charm (at RHIC) is via gluon-gluon interactions, 
and thus the yields at forward rapidity, the deuteron- 
going direction, are sensitive to low-x in the gold nucleus 
where gluon shadowing 0,0] and/or gluon saturation ef- 
fects [5| become important, providing a crucial test of 
these effects. 



There is also significant interest in determining the 
color screening length in the quark-gluon plasma for tem- 
peratures T > 170 MeV, as achieved in relativistic heavy 
ion collisions @ . One proposal for determining this is the 
measurement of several quarkonia states where the dif- 
ferent binding energies (and thus radii) can bracket the 
screening length of interest [7|, |8|. However, this suppres- 
sion of quarkonia must be separated from the aforemen- 
tioned cold nuclear matter effects. Thus precise measure- 
ment of quarkonia suppression in d+Au is needed. 

The PHENIX experiment has previously published 
J/ip results in d+Au collisions at y/s NN = 200 GeV @ 
from data taken in 2003. In this paper we present results 
from d+An collision data taken in 2008, representing an 
increase in yield by a factor of 30-50 over the previous 
results and a reduction in the systematic uncertainties by 
up to a factor of two. Additionally, the p + p reference 
data sets are updated to include higher statistics data 
from 2006 and 2008. 

The PHENIX apparatus is described in detail in [lo| . 
It consists of two sets of spectrometers referred to as the 
central arms, which measure single-particles emitted over 
pscudorapidity (|t^| < 0.35), and the muon arms, measur- 
ing single muons over pseudorapidities (1.2 < \r/\ < 2.4). 
J / if) particles are measured via their dielectron (dimuon) 
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decays at mid (backward and forward) rapidities, and de- 
tailed analysis methods are given in 0, . The d+Au 
data used for this analysis were recorded using selective 
Level- 1 triggers in coincidence with a minimum bias in- 
teraction requirement, which requires one hit in each of 
two beam-beam counters (BBCs) located at positive and 
negative pseudorapidity (3 < \r/\ < 3.9). This minimum 
bias selection covers 88 ± 4% of the total d+Au inelas- 
tic cross section of 2260 mb [Hj]. This can be corrected 
to an unbiased sample, 100% of the total cross section, 
by applying a bias correction factor (c = 0.89) to the 
particle yield measured in any minimum bias event. Ad- 
ditional Level- 1 triggers independently require (1) one 
hit above threshold (600 or 800 MeV) in the Electromag- 
netic Calorimeter with a matching hit in the Ring Imag- 
ing Cerenkov Detector identified as an electron or (2) two 
tracks identified as muon candidates Q. The data sets 
sampled via the Level- 1 triggers represent analyzed inte- 
grated luminosities for the different spectrometers of 54.0 
nb _1 to 69.3 nb _1 . We use p + p reference data for the 
midrapidity dielectrons from (131 ] . For the forward and 
backward rapidity dimuons, we report here new p+p data 
from 2006 and 2008 with a total integrated luminosity of 
5.1 pb~ x . 

The pT-integrated J/ip yield as a function of rapidity 
is calculated via: 



dN 



cN 



j/ij, 



dy NuB^AAy 



(1) 



where Bu is the branching fraction for J/ip — > e + e~ or 
2Vjm is the number of J/ip counts, c is the bias 
correction factor, -/Vmb is the number of sampled mini- 
mum bias events, Ay is the width of the rapidity bin, 
and eA represents the product of the efficiency and ac- 
ceptance corrections, including the Level-1 trigger effi- 
ciency. The number of J/ip particles is determined us- 
ing the invariant mass distribution of unlike-sign lepton 
pairs. Approximately 38000, 8900, and 42000 J/ip counts 
are measured at backward, mid, and forward rapidity, re- 
spectively. Figure [1^, shows the J / ip yields in p+p and 
d+Au unbiased collisions. 

We quantify the cold nuclear matter effects by calcu- 
lating the nuclear modification factor RdAu as given by: 



RdAu{i) 



dN d+Au (i) 

dy 

dNP+P 



<JVcoii(i)> 
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where i is the index of the centrality bin and (N co \\(i)) is 
the average number of nucleon-nucleon collisions and is 
determined from the total energy deposited in the BBC 
located at negative rapidity. For a given centrality bin 
(N co ii(i)) is derived using a Glauber calculation coupled 
to a simulation of the BBC response (as described in 
Q). The centrality bins used in this analysis are char- 
acterized as follows: central (7V co n(0-20%)) = 15.1 ± 
1.0, (iV co ii(20-40%)) = 10.3 ± 0.7, (N coll (40-60%)) = 
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FIG. 1: (color online) (a) J/ip B u dN/dy in p+p and d+Au 
collisions as a function of rapidity. The d+Au yields are di- 
vided by the average number of nucleon-nucleon collisions 
(Af co n(0-100%)} = 7.6. (error lines) point-to-point uncorre- 
lated uncertainties; (boxes) point-to-point correlated uncer- 
tainties; (text) global normalization scale uncertainties, (b) 
RdAu nuclear modification factors for unbiased collisions. 



6.6 ± 0.6, (iVcoii (60-88%)) = 3.2 ± 0.2 and unbiased 
(NcoU (0-100%)) = 7.6 ± 0.3. Figure Hb shows the R dAu 
corresponding to unbiased collisions. Figure^ (b) shows 
RdAu corresponding to d+Au centralities of 60-88% (0- 
20%). Note that more central collisions correspond to 
cases where the nucleons in the deuteron strike closer to 
the middle of the gold nucleus, and thus nuclear effects 
are expected to be enhanced (which is seen in the data) . 

The peripheral RdAu favors some suppression at all 
rapidities, though this result is tempered by the current 
systematics of approximately ±15%. The central RdAu 
indicates a much larger suppression for J/ip at forward 
rapidity. 

We also calculate the ratio Rcp as the nuclear mod- 
ification between central and peripheral c?+Au collision 
classes of events: 



dN 



d+Au 



(0-20%) 



Rci 



dy 



/<JVcoii(0-20%)) 



dJvd+Au (6 - 88%) /(jVeoii(60 _ 88%)) 



(3) 



Figure [2J; shows the Rqp ratio for the most central 
category relative to the peripheral 60-88% category as a 
function of rapidity. The quantity Rcp has the advan- 
tage that many of the systematic uncertainties cancel in 
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the ratio. One observes a dramatic suppression of for- 
ward rapidity yields for central d+Au events compared 
to peripheral events. At backward rapidity, there is little 
to no modification seen. 
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FIG. 2: (color online) Nuclear suppression factors (RdAu) for 
(a) peripheral, (b) central, and (c) Rcp as a function of ra- 
pidity. 

We confront two specific calculations available in the 
literature with our data. The first class of calculations 
that are often employed include nuclear-modified PDFs 
and an effective <7t, r . Here we utilize the EPS09 nuclear 



modified PDF set [14[ and a o\, r = 4 mb (chosen to match 
the unbiased backward rapidity RdAu data) [l5[ • We find 
a reasonable agreement within uncertainties with the un- 
biased RdAu data shown in Fig. [TJd. We also show as red 
dashed lines the differences within the EPS09 nPDFs for 
a single parameter change that gives the largest varia- 



tion 14]. However, if one plots the same calculations for 
Rcp as shown in Fig. [2];, one sees reasonable agreement 
at backward and midrapidities, but a significant under- 
prediction of the suppression for forward rapidity J/ ip in 
the most central events. 

A second class of calculations incorporates gluon sat- 
uration effects at small- X 0, [l6| , and is compared with 
experimental data in Figs. [T][2] A modest J/t/j enhance- 
ment is predicted at midrapidity due to double-gluon ex- 
change processes (not seen in the data) and a substantial 
J I ip suppression at forward rapidity and in more central 



d+Au events due to saturation effects (in agreement with 
the data). 

In order to further explore the centrality dependence 
of the nuclear effects we categorize each d+Au central- 
ity class in terms of the distribution of transverse radial 
positions (ry) of the nucleon-nucleon collisions relative 
to the center of the gold nucleus. The tt distributions 
for the four centrality categories are shown in Fig. 
We expect that the nuclear effects are dependent on the 
density weighted longitudinal thickness through the gold 

nucleus (a(tt) = J d z p( z , r T)^j , where po is the den- 
sity in the middle of the nucleus. This quantity is also 
shown in Fig. [5^, as a function of tt- 

We now posit three different functional dependencies 
of the nuclear modification on K(tt)- 



Exponential : M(r T ) = e ~ aA ( rT ) 
Linear : M(r T ) = 1.0 - aA(r T ) 
Quadratic : M(r T ) = 1.0 - aA(r T ) 2 , 



(4) 
(5) 
(6) 



where a is a parameter depending on the average level 
of modification. The EPS09 nPDF based calculation, 
shown in Figs. [I][2] assumes the linear relation 15|, |l7[ in 
Eq. [5] in order to make centrality-dependent predictions. 
In contrast, contributions from a break up of the cc via 
a <Tbr follow the exponential relation in Eq. |4j 

Figure (3Jd shows the nuclear modification Rcp in the 
most central bin versus the (unbiased) average modifi- 
cation RdAu- This particular set of quantities is cho- 
sen because for each of the three geometric dependencies 
(Eqs.[3H6]), a given value of the parameter a results in a 
unique point on the plot and varying the parameter a re- 
sults in a unique locus of points on which any suppression 
with that geometric dependence must lie. 

The experimental data is also plotted in Fig. [3)3 for the 
same quantities. The ellipses represent a one standard 
deviation contour for the systematic uncertainties, which 
are largely uncorrelated between the unbiased i?dAu and 
Rcp- There is a substantial deviation between the ex- 
ponential and linear cases and the experimental data at 
forward rapidity, while at mid and backward rapidities 
the data cannot discriminate between the cases. Thus, 
the forward rapidity data suggest that the dependence 
on A(t't) is non-linear and closer to quadratic. If the 
dominant physics leading to the modification is different 
at different rapidities, it is possible for example that the 
modification at backward rapidities is linear while at for- 
ward rapidities it is not. This is reinforced by the EPS09 
plus Ohr calculation where regardless of the variation of 
the nPDF or <7„ r one cannot simultaneously describe the 
full centrality dependence of the data as seen in Fig. [21 

Other non-linear density effects (e.g., quadratic) for 
the geometric dependence [3] and for break up of the 
cc after production [l9|, [2(| have been proposed. An al- 
ternative explanation is that initial-state parton energy 
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FIG. 3: (color online) (a) Normalized to unity at the maxi- 
mum bin are (solid curves) transverse radial tt distributions 
in the gold nucleus for four d+Au centrality selections and 
(dashed curve) density weighted longitudinal thickness as a 
function of Tt (A(t-t)). (b) (points) Unbiased RdAu versus 
Rcp for the experimental data and (curves) constraint lines 
for three geometric dependencies of the nuclear modification. 



loss results in a backward shift of the J /tp rapidity dis- 
tribution 



21| . It has been observed [22[ that the nuclear 



modification as a function of center-of-mass rapidity is 
similar to that observed at lower energies [l[ with a steep 
increase in suppression at forward rapidities, as predicted 
for initial-state parton energy loss. 

In summary, we have presented precision data on J/ip 
yields in d+Au and p+p collisions at y/s~^ = 200 GeV 
over a broad range in rapidity and rf+Au centrality. Nu- 
clear modification factors at forward rapidity as a func- 
tion of centrality cannot be reconciled with a picture of 
cold nuclear matter effects (nPDFs and a <7{, r ) when an 
exponential or linear dependence on the nuclear thick- 
ness is employed. Effects of gluon saturation may play 
an important role in understanding the forward rapid- 
ity modifications, though other explanations involving 



initial-state parton energy loss need further investigation. 
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